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optical system, generating an initial estimate of a phase of the 
optical system, iteratively performing a phase retrieval on the 
intensity data to generate a phase estimate using an initial 
diversity function corresponding to the intensity data, gener- 
ating a phase map from the phase retrieval phase estimate, 
decomposing the phase map to generate a decomposition 
vector, generating an updated diversity function by combin- 
ing the initial diversity function with the decomposition vec- 
tor, generating an updated estimate of the phase of the optical 
system by removing the initial diversity function from the 
phase map. The method may further include repeating the 
process beginning with iteratively performing a phase 
retrieval on the intensity data using the updated estimate of 
the phase of the optical system in place of the initial estimate 
of the phase of the optical system, and using the updated 
diversity function in place of the initial diversity function, 
until a predetermined convergence is achieved. 
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HYBRID DIVERSITY METHOD UTILIZING 
ADAPTIVE DIVERSITY FUNCTION FOR 
RECOVERING UNKNOWN ABERRATIONS IN 
AN OPTICAL SYSTEM 

5 

CROSS REFERENCE TO RELATED 
APPLICATIONS 

The benefit of priority is claimed under 35 U.S.C. 1 1 9(e) of 
commonly assigned U.S. Provisional Patent Application No. to 
60/7 13,6 1 6, filed Aug. 3 1 , 2005, entitled “Wavefront Sensing 
Devices and Methods,” which is incorporated by reference 
herein in its entirety. 

STATEMENT OF GOVERNMENT INTEREST 15 

The invention described herein was made by an employee 
of the United States Government and may be manufactured 
and used by or for the Government of the United States of 
America for governmental purposes without the payment of 20 
any royalties thereon or therefor. 

FIELD OF THE INVENTION 

The present invention relates generally to optical systems 25 
and, in particular, to the recovery of aberrations in optical 
systems. 

BACKGROUND OF THE INVENTION 

30 

Phase retrieval is an image-based wavefront sensing 
method that utilizes point-source images or other known 
objects to recover optical phase information. 

Conventional image-based phase retrieval techniques may 
be classified into two general categories: iterative-transform 35 
methods and parametric methods. Modifications to the origi- 
nal iterative-transform approach have been based on the intro- 
duction of a defocus diversity function or on the input-output 
method. Various implementations of the focus-diverse itera- 
tive-transform method have appeared that deviate slightly by 40 
utilizing a single wavelength or by varying the placement and 
number of defocused image planes. Modifications to the para- 
metric approach include minimizing alternative merit func- 
tions as well as implementing a variety of nonlinear optimi- 
zation methods such as Levenburg-Marquardt, simplex, and 45 
quasi-Newton techniques. The concept behind an optical 
diversity function is to modulate a point source image in a 
controlled way. The purpose of a diversity function in wave- 
front sensing is to “spread out” the optical response to expose 
aberration information. In principle, any known aberration 50 
can serve as a diversity function, but defocus is often the 
simplest to implement and exhibits no angular dependence as 
a function of the pupil coordinates. Since defocus has only 
angular dependence, no one part of the point spread function 
is emphasized more than any other. 55 

The iterative-transform method can result in phase 
retrieval recoveries with higher spatial frequency content than 
is practical with parametric methods that solve directly for 
aberration coefficients of a given basis set, such as, for 
example, a Zernike basis set. From a different viewpoint, the 60 
method operation may be analyzed as consecutive projections 
onto convex sets. A disadvantage of the iterative-transform 
approach is that recovery results are prone to spurious phase 
values and branch points, as well as wrapped phase values for 
aberration contributions larger than one wave. As a result, it is 65 
well known that the iterative-transform method, while espe- 
cially suited for high-spatial frequency phase recoveries, has 
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a limited dynamic range when compared with the parameter- 
based approach, in which multi-wave recoveries are possible 
but are generally of lower spatial frequency. Multi-wave 
phase recoveries are possible, but generally require an addi- 
tional phase-unwrapping post-processing step. 

One particular approach is the Misell-Gerchberg-Saxton 
(MGS) method. The MGS approach is well suited for high 
spatial frequency phase map recoveries, but is not well suited 
for multi-wave recoveries, due to the well-known 2 jt phase 
wrapping problems inherent in the iterative-transform 
approach. Nevertheless, the MGS method approach can be 
utilized for multi -wave phase recoveries but requires an addi- 
tional phase-unwrapping post-processing step. 

For the reasons stated above, and for other reasons stated 
below, which will become apparent to those skilled in the art 
upon reading and understanding the present specification, 
there is a need in the art for alternative approaches for recov- 
ering aberrations in optical systems. 

SUMMARY OF EXEMPLARY EMBODIMENTS 

Particular embodiments of the method of the present inven- 
tion are presented. It should be understood that the invention, 
in its broadest sense, could be practiced without having one or 
more features of these embodiments. It should also be under- 
stood that these embodiments are merely exemplary. 

In one embodiment, the invention provides a method of 
recovering unknown aberrations in an optical system, com- 
prising collecting intensity data produced by the optical sys- 
tem, generating an initial estimate of a phase of the optical 
system, iteratively performing a phase retrieval on the inten- 
sity data to generate a phase estimate using an initial diversity 
function corresponding to the intensity data, generating a 
phase map from the phase retrieval phase estimate, decom- 
posing the phase map to generate a decomposition vector, 
generating an updated diversity function by combining the 
initial diversity function with the decomposition vector, gen- 
erating an updated estimate of the phase of the optical system 
by removing the initial diversity function from the phase map . 
The method may further comprise repeating the process 
beginning with iteratively performing a phase retrieval on the 
intensity data using the updated estimate of the phase of the 
optical system in place of the initial estimate of the phase of 
the optical system, and using the updated diversity function in 
place of the initial diversity function, until a predetermined 
convergence is achieved. 

In another embodiment, the invention provides a method of 
recovering unknown aberrations in an optical system, com- 
prising introducing one or more known aberrations into the 
optical system, for each known aberration, collecting inten- 
sity data produced by the optical system, generating an initial 
estimate of a phase of the optical system, iteratively perform- 
ing a phase retrieval on each set of the intensity data to 
generate a phase estimate using an initial diversity function 
corresponding to the sets of the intensity data, combining 
each of the phase retrieval phase estimates to generate a phase 
map, decomposing the phase map to generate a decomposi- 
tion vector, generating an updated diversity function by com- 
bining the initial diversity function with the decomposition 
vector, generating an updated estimate of the phase of the 
optical system by removing the initial diversity function from 
the phase map. The method may further comprise repeating 
the process beginning with iteratively performing a phase 
retrieval on each set of the intensity data using the updated 
estimate of the phase of the optical system in place of the 
initial estimate of the phase of the optical system, and using 
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the updated diversity function in place of the initial diversity 
function, until a predetermined convergence is achieved. 

Aside from the structural and procedural arrangements set 
forth above, the invention could include a number of other 
arrangements, such as those explained hereinafter. It is to be 
understood that both the foregoing description and the fol- 
lowing description are exemplary only. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a schematic and photograph of the testbed tele- 
scope optical design and primary mirror used in validating 
embodiments of the invention; 

FIG. 2 is a visualization of phase estimation in accordance 
with an embodiment of the invention; 

FIG. 3 is a visualization of intensity data used for the phase 
estimation of FIG. 2; 

FIG. 4 is a block flow diagram of a method of recovering 
aberrations in an optical system in accordance with an 
embodiment of the invention; 

FIGS. 5a -5c are visualizations demonstrating filtering arti- 
ficial phase artifacts from phase retrieval results for use with 
embodiments of the invention; 

FIG. 6 is a visualization of one example of an orthogonal 
aberration basis set for use in reporting phase retrieval results 
and tracking aberration terms for use with embodiments of 
the invention; 

FIG. 7 is a diagram of the pupil geometry of the testbed 
telescope used in validating embodiments of the invention; 

FIG. 8 is a visualization of orthogonal and non-orthogonal 
basis functions on segment 9 of the pupil geometry of FIG. 7 
for use with embodiments of the invention; 

FIG. 9 a is a visualization of diversity defocus images asso- 
ciated with individual testbed telescope primary mirror seg- 
ments generated from a single source in accordance with an 
embodiment of the invention; 

FIG. 96 is a visualization of a phase retrieval result from the 
diversity defocus image of segment 6 of FIG. 9 a in accor- 
dance with an embodiment of the invention; 

FIG. 9c is a visualization of a wavefront residual result 
resulting from the subtraction of the diversity defocus image 
of segment 6 of FIG. 9 a from the phase retrieval result of FIG. 
9b in accordance with an embodiment of the invention; 

FIG. 19a is a visualization of intensity data of the diversity 
defocus image of segment 6 of FIG. 9 a in accordance with an 
embodiment of the invention; 

FIG. 106 is a visualization of intensity data generated from 
the phase retrieval result of FIG. 9b in accordance with an 
embodiment of the invention; 

FIGS. lOolOd are visualizations of a comparison of inten- 
sity data of cross-sections of the actual data of FIG. 10a to the 
estimated data of FIG. 106 in accordance with an embodi- 
ment of the invention; 

FIG. 11a is a visualization of phase retrieval results from 
the diversity defocus images of all segments of FIG. 9a in 
accordance with an embodiment of the invention; 

FIG. 116 is a visualization of wavefront residual results 
resulting from the subtraction of the diversity defocus image 
of each segment of FIG. 9a from its corresponding phase 
retrieval result of FIG. 11a in accordance with an embodi- 
ment of the invention; and 

FIG. 12 is a visualization of the calculated point-spread 
function upon correction for aberrations recovered in FIG. 
11a in accordance with an embodiment of the invention. 
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DETAILED DESCRIPTION OF EXEMPLARY 
EMBODIMENTS 

In the following detailed description of embodiments of the 
5 invention, reference is made to the accompanying drawings 
that form a part hereof, and in which are shown by way of 
illustration specific embodiments that may be practiced. 
These embodiments are described in sufficient detail to 
enable those skilled in the art to practice disclosed subject 
l o matter, and it is to be understood that other embodiments may 
be utilized and that process changes may be made without 
departing from the scope of the claimed subject matter. The 
following detailed description is, therefore, not to be taken in 
a limiting sense, and the scope of the claimed subject matter 
15 is defined only by the appended claims and equivalents 
thereof. 

Performance of optical systems, such as those found in 
cameras, telescopes, microscopes, and other systems for the 
collection of images, can often be improved by actively cor- 
20 reefing for aberrations in the optical system. Wavefront sens- 
ing is a technique to obtain phase information from a wave- 
front to recover such aberrations. 

In one application, the method of the present invention may 
be utilized on a space-based telescope, such as NASA’s 
25 James Webb Space Telescope (JWST), for example. On 
JWST, the primary mirror is a 6.5 meter segmented mirror 
having 1 8 segments in a hexagonal array. The method accord- 
ing to an embodiment of the invention may be used with an 
image-based wavefront sensing and control system in such a 
30 telescope to align the mirror segments, minimize the effects 
of figure error, and position the secondary mirror of the three- 
mirror anastigmat design. The wavefront sensing method 
specified for JWST is image-based in the sense that point- 
source stellar images are collected and processed to recover 
35 optical phase information. In addition to wavefront sensing 
and control applications with active and adaptive optical sys- 
tems, embodiments of the invention may be used for optical 
system verification and testing, such as, for example, quanti- 
tative measurement of optical system performance. The opti- 
40 cal systems in those applications may comprise at least one of 
a lens, a mirror, an imaging system, a camera, a detector, and 
a laser. The optical systems may also comprise compound 
systems, including combinations of those elements. 

In one example, the method of the present invention has 
45 been used in an image-based wavefront sensing and control 
system in a meter-class laboratory -based testbed telescope. 
Several exemplary embodiments of the method will be 
described with reference to the testbed telescope. 

FIG. 1 shows a schematic and a photograph of the testbed 
50 telescope optical design and primary mirror. While the 
embodiments will be described with reference to their use 
with the testbed telescope, it is noted that the methods 
described herein are applicable to a variety of other optical 
systems, including monolithic mirror systems, compound 
55 lens systems, and other optical systems adapted to generate a 
wavefront in response to a source. 

The various embodiments described herein will be referred 
to as a Hybrid Diversity Method (HDM). According to some 
embodiments, the method produces an unwrapped phase 
60 result, and may be based on a more direct computational 
procedure, thus eliminating additional post-processing steps 
of traditional phase retrieval processes. Further, some 
embodiments may facilitate high-speed multi-processor 
implementations for reduced latency. 

65 The HDM may utilize an iterative-transform phase 
retrieval approach, without focusing on a convergence of the 
iterative-transform phase retrieval or utilizing a phase-un- 
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wrapping post-processing step. Instead, embodiments of the 
HDM utilize a floating, i.e., adaptive, diversity function to 
incorporate feedback during phase recovery to guide the 
phase retrieval process. Further, some embodiments of the 
method may resolve issues associated with conventional 
iterative-transform phase retrieval, including the convergence 
criterion, and unbiased phase estimation, as detailed below. 

Iterative-transform phase retrieval is pixel-based, i.e., it 
operates by iteratively enforcing known constraints between 
conjugate Fourier domains to recover point-by-point phase 
maps. This approach has been termed the “error-reduction” 
method in the literature. 

Various embodiments of the invention may use a variety of 
iterative-transform phase retrieval processes, and are not lim- 
ited to the types discussed above with reference to conven- 
tional methods. 

Embodiments of the invention may overcome limitations 
in the dynamic range of the conventional iterative-transform 
approach because no phase wrapping occurs when a multi- 
wave aberration is incorporated as a known diversity func- 
tion. That is, for focus-diverse phase retrieval according to 
embodiments of the invention, the diversity defocus term is 
often orders of magnitude larger than other unknown aberra- 
tions in the system. Therefore, in those embodiments, if the 
unknown aberrations can be converted to a known diversity 
function during the recovery process, then point-by-point 
phase recoveries of multi -wave aberrations may be realized 
without a phase-unwrapping post-processing step. In some 
embodiments, this is achieved by transferring the resulting 
sensed aberration into a continually updated, i.e., adaptive, 
diversity function during the recovery process. In other 
words, in those embodiments, the method incorporates feed- 
back. The final aberration estimate may comprise the diver- 
sity function itself, after subtraction of any other known diver- 
sity components. As a result, the method may produce a null, 
i.e., a flat, as the aberration content is continually fed into the 
diversity function during subsequent iterations of the method. 

FIG. 2 is a visualization of phase estimation in accordance 
with an embodiment of the invention, illustrating an example 
of the high dynamic range phase estimation that is possible 
with the HDM (1X=633 nm). FIG. 3 is a visualization of 
intensity data used as input for the phase estimation of FIG. 2. 
As depicted in FIG. 2, multi-wave amplitude aberrations, i.e., 
astigmatism of 7.1565 waves and coma of 1.1059 waves, are 
recoverable using embodiments of the invention, producing a 
null with a root mean square (RMS) wavefront error of a small 
fraction of its wavelength. A direct interferometer compari- 
son with an interferogram calculated from the HDM results in 
this example provided agreement to 17 nm RMS at a wave- 
length of 633 nm, or less than 0.5% RMS recovery error. 
Although a point-source was used to generate the various 
examples shown herein, other sources are suitable for use 
with embodiments of the invention, such as extended known 
sources and extended unknown sources. As used herein, 
“extended source” means any multi -dimensional source, gen- 
erally not approximating a point- source, whether it is of 
known or unknown dimensions. 

As discussed above, this nulling behavior of some embodi- 
ments of the method in turn may provide a natural conver- 
gence criterion, since the measure of convergence may now 
be based on the RMS wavefront error of the null (similar to 
specifying the figure quality of an optical reference flat). 
Being adaptive in nature, some embodiments of the method 
may reduce the estimation bias that stems from imperfect 
prior knowledge of the diversity fimction. That is, diversity 
defocus is not perfectly known at the outset, but a refined and 
corrected estimate may be obtained during each iteration of 
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the method. In one embodiment, the initial diversity function 
and/or phase estimate can be presumed to be some predeter- 
mined value. In one example, the predetermined value may be 
zero. Other non-zero values may also be used. In some 
5 embodiments, the initial phase estimate comprises the output 
of another method, such as an estimation method, for 
example. 

FIG. 4 is a block flow diagram of a method for recovering 
aberrations in an optical system in accordance with an 
to embodiment of the invention. An initial estimate of the phase 
is used to start the method, which is designed to process one 
or more intensity images simultaneously that differ in diver- 
sity defocus levels. These defocus levels may be chosen 
according to well known guidelines. However, such selection 
15 is not critical to performing the embodiments of the method 
described herein, as the embodiments may be performed with 
arbitrarily chosen levels of aberration. And while improved 
accuracy may be facilitated by using multiple images, the 
embodiments are suitable for use with only one set of inten- 
20 sity data, i.e., one wavefront image. 

The embodiment of the invention described above utilizes 
focus diversity, in which a known defocus aberration is intro- 
duced into an optical system. Other embodiments may utilize 
other known aberrations, including, for example, one or more 
25 of piston, tip, tilt, defocus, astigmatism, coma, spherical, and 
trefoil. 

As shown in FIG. 4, the individual images are then pro- 
cessed to obtain phase information. For example, the images 
may be processed using a nonlinear optimization procedure 
30 or a modified MGS iterative transform method. Some 
embodiments of the method are model -based in the sense that 
method parameters, such as defocus, 17#, wavelength, and 
detector pixel size, and data parameters, such as image bias 
and flux, are estimated directly from the image data as a 
35 pre-processing step prior to phase retrieval. This process may 
be handled, for example, using a nonlinear least squares or an 
optimization method such as the Levenburg-Marquardt 
method. In practice, the values of these parameters may affect 
how the Fourier transforms are implemented in the iterative- 
40 transform method. 

As shown in FIG. 4, the output of the iterative-transform 
procedure at this stage is a phase estimate for each image of 
the defocus diversity data set. Any phase wrapping issues of 
45 the iterative-transform procedure may be ignored. These indi - 
vidual phase estimates are then combined, such as by using a 

weighted average (denoted by the ® symbol), to form the 
phase estimate, ty estt „ ew , as illustrated in FIG. 4 and labeled as 
the “inner-loop.” It is well known that the iterative-transform 
50 method stagnates quickly after a relatively small number of 
iterations. Therefore, the overall method efficiency is 
improved from a computational standpoint by restricting the 
inner-loop iterations to a small number. In the implementa- 
tion considered here, this process is repeated for three itera- 
55 tions to take advantage of the rapid initial gain in convergence 
that is characteristic of the iterative-transform method, while 
avoiding stagnation. Other numbers of iterations may also be 
used. Because of the feedback of the updated diversity func- 
tion into subsequent iterative-transform phase retrieval cal- 
60 culations, the need to force convergence of any one iterative- 
transform procedure through repeated iterations may be 
obviated. 

After exiting from the inner loop, a phase map is produced 
that is partially consistent with each image used in the data 
65 set. This phase map is then decomposed into a basis set for 
representation of wavefront data. A basis set may be selected 
that is suitable for the aperture shape. The basis set may 
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comprise orthogonal or non-orthogonal functions. In one 
example, the phase map may be decomposed into a Zemike 
basis set for a given sampling and aperture shape. In another 
example, a fully orthogonal basis set may be used. The vector 
of decomposition coefficients may be written as: 

a diversity&O rthogonal Fit^J. Eq. (1) 

This process is shown in FIG. 4 where a =(a l5 a 2 , . . . , a nz ) 
is the decomposition vector and nz =1 / 2 (z order + 1 )(z order +2) is 
the number of terms. These numerical values are then trans- 
ferred to the previous diversity function using the update: 

^ diversity,!? ~ ® diversity,??- ® diversityjit' Eq. (2) 

This process may incorporate higher-order basis sets in the 
aberration off-loading process. Alternatively, the order four 
set (nz=l 5) may be utilized to capture the controllable modes 
of an optical system. 

In the example implementation global tip and tilt corre- 
spond to the coefficients a 2 and a 3 , which are used to update 
the image registration of individual diversity images. An 
updated diversity data set may then be produced for re-pro- 
cessing by the inner-loop via return to the outer-loop. In this 
embodiment, image registration may be adaptive and may be 
based on the amount of tip and tilt returned in the outer-loop 
iterations. Image-registration errors which may be present 
due to the image centering pre-processing step may be thus 
averaged out as the global tip and tilt are minimized by the 
method. Similarly, in some embodiments, imperfect knowl- 
edge of either the image f/# or diversity defocus value for 
individual images can lead to the recovery of an overall global 
defocus error. However, diversity data may optionally not be 
updated in performance of the methods. 

This possibility of recovering a global defocus error is due 
to an ambiguity between f/# and defocus 6 z, which are related 
in the paraxial region by 

6z=8Af # 2 a rfc/ , Eq. (3) 

where a de As the defocus term in quadratic curvature units, i.e. 
^dep^dej P 2 where p is the normalized radial pupil coordinate. 
In the Zernike basis set the defocus coefficient, a 5 , is related 
to the quadratic curvature coefficient by a factor of two: 

a 5=aV 2 - 

In the illustrated embodiment of the method, previous 
diversity defocus values are adaptively updated according to 
the recovery of global defocus in the phase estimate, § estnew - 
In turn, as the amount of diversity defocus varies in each 
image, aberration recovery varies with differing amounts of 
defocus and, thus, feedback is incorporated into the recovery 
process. As shown, this process is continued until the global 
defocus error is driven to some predetermined convergence 
criteria, and may be driven substantially to zero. 

In some embodiments, each aberration contribution may 
far exceed one wave after exiting the inner-loop portion of the 
method and, as is well known, the classical iterative-trans- 
form method cannot, by itself, recover multi-wave aberra- 
tions . Accordingly, the method may produce a wrapped phase 
map, ^ in the absence of some form of aberration 
“off-loading.” However, an aberration fitting procedure may 
still be applied in Eq. (1) to this wrapped phase data to return 
at least some portion of the multi-wave aberration to the 

variable a diversity ^ v These numerical values may then be 
transferred to the diversity function and become part of the 
known phase values . In this way, a multi -wave aberration may 
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be recovered incrementally by successively applying the 
aberration fitting procedure to intermediate wrapped phase 
maps. 

As shown in FIG. 4, during recovery, as more and more of 
5 the aberration is transferred to the diversity function follow- 
ing successive iterations of the outer loop, the estimated 
phase, <fy estf „ ew , no longer wraps since the multi-wave (un- 
known) aberration values become incorporated as part of the 
diversity function. As a result, § est ^ new begins to resemble a 
to reference flat, i.e., a null (see FIG. 2), as the aberration content 
is transferred to the diversity function. The process of aber- 
ration transfer is illustrated in FIG. 4 with a calculated inter- 
ferogram from both the adaptive diversity function, 

15 a diversity fiv and the reference flat, <p estt „ ew , as the method 
nears convergence for a wavefront dominated by astigma- 
tism. 

In some cases when representing obscurations or edges in 
the pupil constraint of the phase retrieval process, the Fourier 
20 transform used in diffraction modeling may alias power into 
the complex amplitude of the image plane point spread func- 
tion (PSF). This artificial power may be problematic when 
using coarsely sampled edges and can contribute to an artifi- 
cial and undesirable “ringing” that appears in the recovered 
25 phase maps. To reduce this phase artifact, the various embodi- 
ments may use a filtering technique based on Gaussian 
smoothing, to soften edges in a coarsely sampled aperture. 
Other techniques also exist, such as, for example, edge 
weighting based on the area of pixel coverage along a seg- 
30 ment edge, and may be utilized with the various embodi- 
ments. FIGS. 5 a -5c are visualizations demonstrating filtering 
artificial phase artifacts from phase retrieval results for use 
with embodiments of the invention. FIG. 5 a shows an 
example of aliasing from a fast Fourier transform (FFT) of a 
35 hexagon mask without filtering. FIG. 5b shows the analytic, 
i.e., ideal, response of the same hexagon mask. FIG. 5c shows 
an example of aliasing from a FFT of a Gaussian- smoothed 
hexagon mask. 

In some cases, the Gaussian- filtered result may compare 
40 more favorably to the analytical result, although some 
residual aliasing effects still appear. One may conclude from 
FIGS. 5a -5c, that a basic trade occurs between aliased power 
and the spatial frequency response of the FFT method. The 
implication for the iterative-transform process is that the mid- 
45 or higher- spatial frequency components may be damped in 
the recovered phase maps. This smoothing or filtering may 
effectively reduce artificial phase artifacts due to aliasing. 

As discussed above, orthogonal or non-orthogonal basis 
sets may be used with embodiments of the method of the 
50 invention. A basis set may be chosen that is best suited to the 
aperture under consideration. For example, the Zemike basis 
functions are orthogonal over a circular aperture. An alterna- 
tive basis set may be developed for more complex shapes, 
such as, for example, the hexagonal apertures used herein for 
55 the example implementation with the testbed telescope mir- 
ror segments. Basis functions for hexagonal apertures, called 
“hexikes,” have been previously derived. The order 4 modes 
(15 terms) are illustrated in FIG. 6 . 

The various embodiments are further adaptable to more 
60 complex geometries. For example, the testbed telescope pupil 
consists of more complicated apertures than the simple hex- 
agonal shapes of the mirrors (see FIG. 7). Out of 18 testbed 
telescope segments, only five are regular hexagons, due to 
vignetting that results from alignment fiducials and second- 
65 ary mirror safety support ears. The remaining segment aper- 
tures deviate significantly from hexagonal shapes. A set of 
orthogonal basis functions may be generated over a sampled 
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irregular aperture starting from a set of “seed” basis func- 
tions, based on well-known Gram-Schmidt orthogonaliza- 
tion. FIG. 8 shows a comparison of how several of the 
orthogonal basis functions compare with the standard 
Zemike basis functions for segment 9 (selected arbitrarily). 
Such methods are applicable to other non-circular apertures. 

FIGS. 9a-9c correspond to unstacked image data of a nega- 
tively defocused image. The hexagonal boundaries, have 
nothing to do with the testbed telescope pupil, but are simply 
an efficient way of extracting segment data from the image 
array with minimum overlap between neighboring images. 

The residual, i.e., null, wavefront that is produced by HDM 
for segment 6, in accordance with one embodiment, as 
detailed above, is also shown in FIG. 9c. Given that this 
residual is essentially what is left after the first 15 Zemike 
basis functions are extracted, i.e., “off-loaded,” to the diver- 
sity array, stress points on the mirror are clearly visible, 
illustrating the hexagonal symmetry of the hexapod mount, in 
addition to a mid-spatial frequency residual that is a product 
of the mirror figuring process. The dominant aberration term 
for segment 6 is astigmatism. The overall peak-valley wave- 
front error (in double pass) is 1 .85 waves at a wavelength of 
Z=6 33 nm. The residual RMS wavefront error is 45 nm, 
which translates to a method convergence stability of 30 
RMS in single-pass. This means that when using an order 4 
aberration basis set (15 terms) for the aberration transfer 
described above, the aberration transfer has converged to the 
Z/30 level. Different optical systems may have different con- 
vergence results. 

Using the phase-retrieval results shown in FIG. 9b, a model 
of the intensity data is constructed and then compared to the 
actual data. This comparison is made in FIGS. lOa-lOd, 
showing full PSF s and x-y cross sections for both the data and 
model. The double-pass phase-retrieval results shown in FIG. 
9b and FIGS. lOa-lOd are typical of the other testbed tele- 
scope segments. A summary of these results is shown in FIG. 
11 a, indicating that the dominant aberration on the segments 
is astigmatism. 

A visualization of wavefront residual results from the sub- 
traction of the diversity defocus image of each segment of 
FIG. 9a, e.g., the first 15 Zernike basis functions, from its 
corresponding phase retrieval result of FIG. 11a. The results 
of FIG. 116 may serve as a diagnostic, illustrating mirror 
stress points induced by the radius of curvature actuation that 
is used to correct segment focus errors in the example imple- 
mentation. Mirror segments requiring the greatest radius of 
curvature control exhibit a slight edge effect along the interior 
edges of these segments. Also visible is a mid-spatial fre- 
quency residual that is a product of the mirror figuring pro- 
cess. As calculated above, the single-pass residual wavefront 
error for FIG. 116 is approximately 20 nm RMS or X/32, 
establishing the magnitude of the residuals due to actuation 
stress and the mirror figuring process. 

Upon correction for aberrations recovered in FIG. 11a, a 
representation of the source producing the wavefront may be 
calculated. FIG. 12 is a visualization of the calculated point- 
spread function corresponding to the diversity focus images 
of FIG. 9a showing good correlation to a laboratory point- 
source as utilized in generating the example diversity focus 
images. 

While the various embodiments have been described with 
reference to a narrowband or monochromatic wavefront, the 
embodiments of the invention are also suitable for use with 
broadband data, i.e., white light. Traditional iterative-trans- 
form phase retrieval approaches have required such narrow- 
band, monochromatic images. However, narrowband filters 
may reduce signal level and increase detector integration 
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times. Because of the adaptive diversity function of the vari- 
ous embodiments, recovery of phase information and aberra- 
tions may be possible without requiring narrowband filtering. 
In addition to permitting the use of broadband filters, which 
5 produce less signal loss and facilitate faster integration times 
when compared to narrowband filters, the various embodi- 
ments may even be used without filtering of the wavefront, 
thus further reducing system complexity. 

The various embodiments described herein include meth- 
10 ods of recovering unknown aberrations in a variety of optical 
systems. Some embodiments use feedback through an adap- 
tive diversity function and may avoid a phase-unwrapping 
post-processing step typical of conventional iterative-trans- 
form phase retrieval methods. In some embodiments, this 
15 may yield a more direct computational procedure and may 
eliminate one or more additional software processes. 

It will be apparent to those skilled in the art that various 
modifications and variations can be made to the structure and 
methodology described herein. Thus, it should be understood 
20 that the invention is not limited to the examples discussed in 
the specification. Rather, the present invention is intended to 
cover modifications and variations. 

What is claimed is: 

25 1 . A method of recovering unknown aberrations in an opti - 

cal system, comprising: 

collecting intensity data produced by the optical system; 

generating an initial estimate of a phase of the optical 
system; 

30 iteratively performing a phase retrieval on the intensity 
data to generate a phase estimate using an initial diver- 
sity function corresponding to the intensity data; 

generating a phase map from the phase retrieval phase 
estimate; 

35 decomposing the phase map to generate a decomposition 
vector; 

generating an updated diversity function by combining the 
initial diversity function with the decomposition vector; 

generating an updated estimate of the phase of the optical 
40 system by removing the initial diversity function from 
the phase map; and 

repeating the process beginning with iteratively perform- 
ing a phase retrieval on the intensity data using the 
updated estimate of the phase of the optical system in 
4 place of the initial estimate of the phase of the optical 

system, and using the updated diversity function in place 
of the initial diversity function, until a predetermined 
convergence is achieved. 

2. The method of claim 1, wherein the intensity data is 
generated in response to a source selected from the group 
consisting of a point source, an extended known source, and 
an extended unknown source. 

3. The method of claim 1, wherein the intensity data com- 
55 prises broadband data. 

4 . The method of claim 1 , further comprising introducing a 
known aberration into the optical system prior to collecting 
intensity data. 

5 . The method of claim 4 , wherein the known aberration 
60 comprises at least one of piston, tip, tilt, defocus, astigma- 
tism, coma, spherical, and trefoil. 

6 . The method of claim 4 , wherein the known aberration is 
at least an order of magnitude greater than a largest expected 
unknown aberration. 

65 7. The method of claim 1, wherein generating an initial 

estimate of a phase of the optical system comprises presum- 
ing the phase to have a predetermined value. 
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8. The method of claim 1, wherein decomposing the phase 
map comprises decomposing the phase map into a basis set 
suitable for a given sampling and aperture shape of the optical 
system. 

9. The method of claim 1, wherein the predetermined con- 
vergence is a residual RMS wavefront error of less than a 
wavelength of the intensity data. 

10. A method of recovering unknown aberrations in an 
optical system, comprising: 

introducing one or more known aberrations into the optical 
system; 

for each known aberration, collecting a set of intensity data 
produced by the optical system; 

generating an initial estimate of a phase of the optical 
system; 

iteratively performing a phase retrieval on each set of the 
intensity data to generate a phase estimate using an 
initial diversity function corresponding to the sets of the 
intensity data; 

combining each of the phase retrieval phase estimates to 
generate a phase map; 

decomposing the phase map to generate a decomposition 
vector; 

generating an updated diversity function by combining the 
initial diversity function with the decomposition vector; 

generating an updated estimate of the phase of the optical 
system by removing the initial diversity function from 
the phase map; and 

repeating the process beginning with iteratively perform- 
ing a phase retrieval on each set of the intensity data 
using the updated estimate of the phase of the optical 
system in place of the initial estimate of the phase of the 
optical system, and using the updated diversity function 
in place of the initial diversity function, until a predeter- 
mined convergence is achieved. 
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11. The method of claim 10, wherein the intensity data is 
generated in response to a source selected from the group 
consisting of a point source, an extended known source, and 
an extended unknown source. 

5 12. The method of claim 10, wherein the intensity data 

comprises broadband data. 

13. The method of claim 10, wherein the one or more 
known aberrations are chosen from piston, tip, tilt, defocus, 
astigmatism, coma, spherical, and trefoil. 

10 14. The method of claim 10, wherein the one or more 

known aberrations comprise defocus. 

1 5 . The method of claim 10, wherein combining each of the 
phase retrieval phase estimates comprises generating a 
weighted average of the individual phase retrieval phase esti- 

15 mates. 

1 6 . The method of claim 10, wherein combining each of the 
phase retrieval phase estimates further comprises filtering the 
phase retrieval phase estimates prior to generating the phase 
map. 

20 17. The method of claim 10, wherein generating an initial 

estimate of a phase of the optical system comprises presum- 
ing the phase to have a predetermined value. 

18. The method of claim 10, wherein decomposing the 
phase map comprises decomposing the phase map into a basis 

25 set suitable for a given sampling and aperture shape of the 
optical system. 

19. The method of claim 10, wherein the predetermined 
convergence is a residual RMS wavefront error of less than a 
wavelength of the intensity data. 

30 20. The method of claim 10, further comprising: 

generating updated intensity data using at least some com- 
ponents of the updated diversity function; and 

comparing the updated intensity data with the collected 
intensity data to validate convergence. 



